Superconducting gap in the presence of bilayer splitting in underdoped 

(Pb,Bi)2Sr2CaCu208+5. 
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The clearly resolved bilayer splitting in ARPES spectra of the underdoped Pb-Bi2212 compound 
rises the question of how the bonding and antibonding sheets of the Fermi surface are gapped in the 
superconducting state. Here we compare the superconducting gaps for both split components and 
show that within the experimental uncertainties they are identical. By tuning the relative intensity 
of the bonding and antibonding bands using different excitation conditions we determine the precise 
k-dependence of the leading edge gap. Significant deviations from the simple cos(fc2;)-cos(fcy) gap 
function for the studied doping level are detected. 

PACS numbers: 74.25.Jb, 74.72.Hs, 79.60.-i, 71.18.+y 



One of the crucial steps in the evolution of our un- 
derstanding of superconductors has been the discovery 
of the energy gap between the ground state and the 
quasiparticle excitations of the system. Photoemission 
is one of those experimental techiques which provided 
direct evidence for the existence of such aj^s^ in the 
high-temperature superconducjiiig cupratesooB. Sub- 
sequent ARPES investigationgaOla clearly demonstrated 
the anisotropic character of the gap which now is con- 
sidered as a hallmark of the high temperature supercon- 
ductors. Detailed studies of the symmetry pppperties of 
the gap suggested a purely rf^wave characteiQ in case of 
the optimallyEl and overdopedEI samples, as well as notice- 
able deviations from the simplest d-wave scenario upon 
underdopingjj. 

More recently, ARPES has made another step for- 
ward - technical improvements of the method al- 
lowed to resolve additional features in the momentum 
and energy distributions of the photoelectrons which 
were not resolved before. Among the new achiev- 
ments is the recent observation of the bilayer splitting 
(BS) in over- and underdoped Bi2Sr2CaGiJ|aJQg-(,5 and 
(Pb,Bi)2Sr2CaCu208+5 bilayer compoundfMj-B The 
presence of the BS implies that some of the conclusions 
based on the earlier ARPES data should be revised. As a 
result of such a revisit it was already established that, for 
example, the famous broadening of the ARPES features 
in momentum^ and energyt3 could be at least partially 
accounted for by the bilayer splitting. Very recent inves- 
tigations, in which differentiation between bonding and 
antibonding components is essential, imply that the more 
active part of the Brillouin zone in a sense of a coupling 
to other degrees of freedom is the antinodal regimunf the 
Fermi surface rather than the (7r,0)-point itselioEj. At 
the same time detection of the BS not only provides an- 
swers but also results in further new questions. One of 
such natural and straightforward questions in the light of 
these findings is 'how does the superconducting gap be- 



have on both Fermi surface sheets?' Could it be that the 
superconductivity in Bi2212 is a sheet|dependeut.one as, 
for instance, is the case in 2H-NbSe2tJ or MgB2llJ? Does 
the strong correlation between the so-called range param- 
eter (reLated to the hopping integral t') and Tcmax, found 
recentljilj for nearly all hole-doped one-layer and bonding 
subband of multilayer HTSC cuprates, really mean that 
the bonding component plays a more important role in 
superconductivity than the other? 

In this Communication we make a start at answering 
these questions by measuring the superconducting gap 
in modulation-free Pb-Bi2212 single crystals. Exploiting 
the relative variation of the photoionisation matrix ele- 
ments for the bilayer split components, we are able to 
independently determine the leading edge gaps (LEG) 
corresponding to the different sheets of the Fermi surface 
in the superconducting state. Our results suggest that 
both gaps are comparable as follows from the analysis 
of the "map of gaps" around the (7r,0)-point and other 
points in the Brillouin zone. Choosing the excitation con- 
ditions in a way which results in the suppression of the 
photocurrent corresponding to the antibonding compo- 
nent we determine the precise k-dependence of the LEG 
for its bonding counterpart for which emission then is 
favorable. The observed anisotropy of the leading edge 
gap for the studied doping level deviates from the conven- 
tional behavior described by the simple cos(fca;)-cos(fcy) 
function. 

The ARPES experiments were carried out using angle- 
multiplexing electron energy analyzers (SCIENTA SES- 
200 and SES-100). Spectra were recorded either using 
hi/ =21.218pe.V photons from a He source, as described 
elsewherc£3o, or using radiation from the U125/1-PGM 
beamlinecj at the BESSY synchrotron radiation facil- 
ity. The total energy resolution ranged from 10 meV 
(FWHM) at hiy = 21 eV to 20 meV at = 50 eV using 
the linearly polarized light as determined for each excita- 
tion energy from the Fermi edge of polycrystalline gold. 
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FIG. 1: Energy distribution map recorded along the az- 
imuthal path (|k| = 1.09 A"^) at T=40K (upper left panel) 
showing the clearly resolved bilayer splitting. EDO's from 
the k-points marked by the white dotted lines on the EDM 
(lower panels). Comparison of EDO's 2 and 3 on an expanded 
energy scale (upper right panel). 



The angular resolution was kept below 0.2° both along 
and perpendicular to the analyzer entrance slit. Data 
shown in Fig. 4 were taken with 0.2° x 0.3° angular resolu- 
tion. All data were collected on two similar underdoped, 
modulation-free single crystals of Pb-Bi2212 {Tc—77K). 
To achieve underdoping the as-grown samples were an- 
nealed in Ar atmosphere at 450° C during 24 hours with 
subsequent fast cooling. 

We start by presenting in Fig. |l| what motivated us to 
initiate the present study. The explicitely seen bilayer 
splitting in the energy distribution map (EDM) recorded 
along the azimuthal path (see dashed line in the inset 
to Fig. 1^) apparently signals the possibility of existence 
of the two gaps related to the bonding and antibonding 
Fermi surfaces. In the lower panels of Fig. l] we show 
three energy distribution curves (EDC's), which corre- 
spond to cuts located at the numbered white dotted lines 
on the EDM. Spectra 1 and 3 are taken at what would 
be vectors if the system was ungapped and, despite 
being very close to one another in k-space, they have 
completely different lineshapes from which a supercon- 
ducting gap value can be derived. EDC number 2 is 
from just in between the gapped sheets of the FS and 
therefore is not suitable for the determination of the su- 
perconducting gap. The upper right panel of Fig. |l| com- 
pares EDC's 2 and 3 in detail to illustrate how erroneous 
the determination of the superconducting gap could be if 
the precise k-location is uncertain. In this case the peaks 
of the EDC's are separated by ~10meV and the leading 
edges by ~6meV. Depending on which method is used 



for the determination of the gap value (see below), the 
error could amount to as much as 50% of the gap value 
reported before for this k-space region. It is also easy to 
predict an error which is the result of the non-resolved 
BS due to an insufficient angular resolution. Summing up 
EDC's from the larger k-area would inevitably result in 
an overestimation of the gap values since both peak and 
leading edge positions of "irrelevant" EDC's are always 
at higher binding energies. 

The necessity, then, to determine the superconduct- 
ing gaps of both Fermi surface sheets is obvious. The 
only question remaining is 'how best to do this?' The 
problem here is that the situation depicted in Fig. |^ is 
rather unique - it is in fact quite unusual to see bonding 
and antibonding bands equally well resolved. To achieve 
such a clarity, EDM shown in Fig. |] was recorded along 
the specific path in k-space: both Fermi surfaces are cut 
at approximately right angles at points which are suf- 
ficiently away from the node (where BS vanishes) and 
antinodes where antibonding saddle point is close to the 
Fermi level and, due to the underdoping, superconduct- 
ing gaps are large hampering the visual resolution of the 
BS. In most cases excitation conditions which define the 
distribution of the photocurrcnt via matrix eleinents are 
such that the bilayer splitting is hardly visibleEStJ. This 
is most likely the reason why the BS remained elusive 
for so long. Every change of a given excitation energy 
and geometry of the experiment leads to a particular, as 
a rule very strong, k-dependence of the matrix elements 
which is, generally speaking, not the same for the bilayer 
split bands. It has been shown recently that the key fac- 
tor defining relative bonding/ai it i b opding intensity is the 
energy of the exciting photonai3ll3'Eil. We exploit these 
observations here and search for suitable excitation con- 
ditions by varying the energy of the tunable synchrotron 
radiation. The most convenient region of the k-space 
for such an exercise (given modulation-free crystals) is 
the vicinity of the (7r,0)-point containing the antinodal 
points of the normal state Fermi surface: here both the 
gap and the bilayer splitting are expected to be largest. 

In Fig. 1^ we present azimuthal momentum distribu- 
tion curves (MDC) for the path indicated in the inset 
by the double-headed arrow. The MDC signal represents 
an integration over 1 degree of polar angle and the first 
20 meV in binding energy. As expected, variations in the 
lineshapes are dramatic. If using linearly polarized radia- 
tion of hi^=21eV or h.v=38eV emissions from the bonding 
and antibonding bands are comparable, for hz^=25eV and 
hj/=50eV the dominance of antibonding component is ob- 
vious. At first sight, it may appear surprising that the 
MDC's taken with hi/==21eV and hi/=21.2eV are so dif- 
ferent. The two important differences here are polariza- 
tion and energy. Firstly, for the former, linearly polarized 
light was used (E||k|| ), whereas the latter involved major- 
ity unpolarized radiation from a helium resonance source. 
Secondly, both theoretical and experimental studies have 
shown that there is a sharra local minimum in the matrix 
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FIG. 2: Azimuthal MDC's showing how relative intensity of 
the bonding and antibonding components could be tuned by 
choosing a suitable excitation energy. Inset illustrates which 
parts of k-space were probed. Vertical dotted line show the 
angular positions of the bonding and antibonding antinodal 
points. 



even a small change in energy could be enough to provoke 
such a large variation of the photoemission signal. 

At this stage we have to decide which energy would be 
the most suitable for the determination of the supercon- 
ducting gaps related to both Fermi surfaces. This choice, 
in turn, depends on how we plan to extract the gap from 
the spectra. There are two widely-known approaches. 
The first involves fitting an EDC with a model spectral 
function multiplied by the Fermi function and corpioluted 
with the energy and angular resolution functionf£3. The 
advantage of this method is that if the fit is successful, 
one gains access to the absolute value of a gap. The 
shortcomings of this method include increased error bars 
due to the fitting procedure itself and, more importantly, 
a lack of knowledge as regards the precise form of the 
spectral function, thus making this approach model de- 
pendent. The second possibility is the so-called leading 
edge gap method, in which the lowest binding energy 
where the ki?-EDC looses half of its maximal intensity 
is determined. This approach is attractive in its preci- 
sion (nowadays better than 1 meV) and simplicity and 
is by far the most widely apphcd in the HTSC ARPES 
literature. The drawback is that the LEG is only a good 
qualitative measure of the gap and does not yield its ab- 
solute value. Now both of these gap methods face a se- 
vere challenge when two separate, gapped features are in 
close proximity. In this context, the EDC fitting proce- 
dure would benefit from having both features of equal 
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FIG. 3: Left panel: map of gaps derived from the EDC's 
taken with 25 eV photons at k-points within the shaded area 
shown in the inset to Fig.^ Both of the gapped Fermi surface 
sheets are seen. Right panel: leading edge gaps as a function 
of the absolute k-value. Zero gap corresponds to the binding 
energy (-7 meV) of the leading edge of the nodal k^- EDC. 



intensity (as, for instance, is the case when using linearly 
polarized 21eV photons, see Fig. H), whereas the LEG de- 
termination would work better when the intensity from 
the antibonding component dominates, as then the weak 
bonding component will not infiuence the position of the 
antibonding leading edge unduly (see also Fig. |^). As 
the aim of this paper is a comparative study of the rela- 
tive sizes and symmetries of the superconducting gaps of 
the two Fermi surface sheets, the LEG approach is the 
method of choice and thus the 25 eV data are the most 
appropriate to address this question. 

Having chosen the excitation energy and the gap de- 
termination method we are left with the known problem 
of which particular EDC has to be used to define the 
LEG (see Fig. |l|). Here we apply a mapping procedure 
allowing to cover a relatively large k-space region and 
introduce thereby a new method of presenting informa- 
tion derived from ARPES spectra - "the map of gaps" - 
which is the momentum (or angular) distribution of the 
binding energies of the leading edges. 

In Fig. ^ we show such a map which is a grey scale 
image of leading edge binding energy. Four traces of 
the gapped Fermi surfaces are easily seen, despite the 
blurring of the bonding components as can be expected 
from the relative intensities seen in the 25 eV panel of 
Fig. ||. As a next step, we have determined the lines 
formed by joining the "minimum gap loci"Ej, which are 
then shown superimposed upon the map of gaps. The 
LEG values along these minimum gap loci are plotted in 
the right panel of Fig. ^ Thus, as a first major result, 
these curves unambiguously demonstrate that the super- 
conducting energy gaps on both Fermi surface sheets are 
identical, indicating no severe sheet dependence of super- 
conductivity in this case. 

Now, being certain that the gaps are identical, we can 
take a further step in trying to obtain information about 
their symmetry. To increase the accuracy of the mea- 
surements one can reduce now the infiuence of one of 
the bands assuming that both gaps have also the same 
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k-dependence in the other part of the BZ, i.e. mono- 
tonically decrease when approaching the nodal direction. 
Looking again at Fig. || one sees that possible candi- 
dates are the photon energies of 21.2 and 50 eV at which 
emission from the antibonding or bonding bands, respec- 
tively, is significantly suppressed. As mentioned before, 
every experiment using linearly polarized radiation im- 
plies particular symmetry conditions which uniquely de- 
fine the photocurrent distribution. In cuprates, due to 
the symmetry of the states forming the near-i?^;- elec- 
tronic structure, it is known that favorable conditions 
for the emission along the F — (tt, tt) direction are such 
that the vector of polarization should be perpendicular 
to this direction whereas for F — (tt, 0) direction it has 
to be parallelHES. In our experimental geometry, when 
using the linearly polarized excitation the polarization 
vector is always parallel to the studied direction in the 
reciprocal space and therefore cannot be equally effective 
for all states at the Fermi level. This implies that to be 
sure not to loose the spectral weight due to symmetry 
reasons, which is especially important when determin- 
ing LEG, it is necessary to probe every k-point with two 
mutually perpendicular orientations of the polarization 
vector which is technically quite difficult. For this reason 
we present here data taken with 21.2 eV radiation from 
the helium lamp which is mostly unpolarized. 

Precise determination of k-dependence of the super- 
conducting gap and thus its symmetry requires precise 
navigation in k-space in order to accurately locate a posi- 
tion on the (former) Fermi surfaceEZI. For this purpose we 
show in Figure ^ a map of intensity integrated within 50 
meV energy interval which facilitates identification of the 
k- vectors. We note here that not only this map but also 
an additional Fermi surface map taken at room temper- 
ature covering much larger area in k-space were needed 
to assure accurate definition of the angle (jX^. Fig. ||b 
shows the greyscale "map of gaps". Visual inspection of 
the map of gaps already clearly points to the anisotropic 
character of the gap. 

One can now use two approaches to construct the 
LEG((/3) function. One method would be to define the 
path in the k-space using, for example, maxima of MDC's 
taken perpendicular to the FS and then plot the value of 
the LEG along this path as a function of ip. We show 
here the results of an alternative method which is based 
on the " minimum gap locus" , but point out that both 
LEG{(p) curves qualitatively agree. Adopting this second 
procedure, we simply replot the same dataset (Fig. ^) in 
other coordinates: LEG value versus Fermi surface angle 
( Fig. |c). This gives the red symbols shown in Fig. 
The curve joining the low-gap extremity of these data 
points represents the true k-dependence of the supercon- 
ducting LEG, and is shown as a dotted line in Fig. 
In this way, these extremal LEG values correspond to 
the points on the "minimum gap locus" dotted line in 
Fig. lb). 

As is clear from the comparison shown in Fig. that 
the obtained k-dependence of the LEG for these UD77K 
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FIG. 4: a) Angular distribution map of the photoemission in- 
tensity integrated within 50 meV of Ef showing the locus of 
kF-points. b) Map of gaps derived from the EDO's taken with 
21.2 eV photons within the similar k-space area. Both maps 
are presented in angular space and therefore give a slightly 
unfamiliar view of the Fermi surface, c) Plot of the leading 
edge gaps as a function of the angle around the Fermi sur- 
face. The red symbols are the data from panel (b), dotted 
lines in both (b) and (c) show the minimum gap locus and 
corresponding LEG values respectively, the solid line in (c) is 
the simple d-wave gap function A{(p) = Amax\cos{2ip)\. 



crystals is quite different from that expected in the case of 
a simple d-wave gap function cos(A:a;)-cos(/cj,). The most 
striking difference is that the gap behaves much flatter 
around the node resulting in a "U"-like shape rather than 
"V"-shape for LEG((p) function. Keeping in mind that 
LEG is only qualitative measure of a real gap it is crucial 
to understand to which extend LEG(<p) reproduces the 
main features of a reaLgap function. According to our 
numerical simulationscJ measured binding energy of the 
nodal ki?-EDC leading edge (-7 meV) corresponds to the 
gapless spectral function thus implying the consistency 
of the presented data with the d-wava^ymmetry of the 
superconducting gap. It is also showrO that due to the 
temperature effects the LEG method tracks the real gap 
with an uncertainty of order of 2 meV (for the experi- 
mental conditions in question) when the corresponding 
absolute values of the real gap are of the same order. 
This last remark lives finite the probabilty that a nodal 
cusp can be present in the real gap function (although 
much more obtuse than is required by the modulus of 
the simplest d-wave function) . 

Observed partial suppression of the superconducting 
gap with respect to the simplest d-wave expectations in 
underdoped sample if-iin qualitative agreement with pre- 
vious ARPES resultgj where authors extracted gap val- 
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ues by modeling the spectra taken at a few k-points using 
a simple BCS spectral function without considering BS. 
Our high resolution data provide much more presice ex- 
perimental information for a further quantitative analysis 
and support the results reported in Ref. ^ in a way that 
as it is clear now, using the 22 eV photons (as in Ref. ||) 
one is sensitive mostly to the bonding componentllj and 
determination of the superconducting gap is not, influ- 
enced by the presence of the BS. Our simulationa23 show 
that the relation between the LEG and the real gap func- 
tion is nontrivial and model-dependent and therefore we 
do not make an attempt here to extract quantitative 
information following authors of Ref. ^ but share their 
opinion that the observed|-behavior could be explained 
invoking higher harmonicsE2l of the gap and is probably 
connected with increased range of the pairing interaction 
in u pdeij doped compounds due to reduced screening (see 
alsoEH). In this context we mention related theoretical 
workEj where momentum-dependent scattering from the 
impurities may also lead to extended gapless regions in 
the gap function centered around the nodes of the pure 
d-wave superconductor. 

As another possible scenario we mention finally that 
a preliminary analysis of ARPES spectra from the same 
samples taken above Tc suggests that the k-dependence 
of the pseudogap is similar to that of the superconduct- 



ing gap presented here. However, at this stage, it is too 
early to draw conclusions from this as to what extent the 
normal state and superconducting state gaps may share 
the same microscopic origin. 

In conclusion, we have demonstrated that in under- 
doped, modulation-free Pb-Bi2212 each of the bilayer 
split Fermi surface sheets supports a gap in the super- 
conducting state. We compared the gap values in the 
antinodal region where both are expected to be maximal 
and found that the gap values on both sheets are iden- 
tical within the experimental error bars. By varying the 
experimental excitation conditions we could tunc the rel- 
ative intensity from the bonding/antibonding states and 
thus, via suppression of the emission from the antibond- 
ing band, we were able to precisely define the symmetry 
of the leading edge superconducting gap for the bond- 
ing Fermi surface sheet. Significant deviations from the 
simplest cos(fc3;)-cos(/cj^) d-wave behavior were observed, 
with a strong flattening out of the gap function around 
the Brillouin zone diagonal. 
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